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ABSTRACT 

Hypervelocity stars (HVS) traverse the Galaxy from the central black hole to 
the outer halo. We show that the Galactic potential within 200 pc acts as a high 
pass filter preventing low velocity HVS from reaching the halo. To trace the orbits 
of HVS throughout the Galaxy, we construct two forms of the potential which 
reasonably represent the observations in the range 5-10 5 pc, a simple spherically 
symmetric model and a bulge- disk-halo model. We use the Hills mechanism 
(disruption of binaries by the tidal field of the central black hole) to inject HVS 
into the Galaxy and to compute the observable spatial and velocity distributions 
of HVS with masses in the range 0.6-4 M . These distributions reflect the mass 
function in the Galactic Center, properties of binaries in the Galactic Center, 
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and aspects of stellar evolution and the injection mechanism. For 0.6-4 M 
main sequence stars, the fraction of unbound HVS and the asymmetry of the 
velocity distribution for their bound counterparts increases with stellar mass. 
The density profiles for unbound HVS decline with distance from the Galactic 
Center approximately as r~ 2 (but are steeper for the most massive stars which 
evolve off the main sequence during their travel time from the Galactic Center); 
the density profiles for the bound ejecta decline with distance approximately as 
r~ 3 . In a survey with limiting magnitude V < 23, the detectability of HVS 
(unbound or bound) increases with stellar mass. 

Subject headings: Galaxy: center - Galaxy: halo - Galaxy: structure - (stars:) 
binaries: general - stars: kinematics - stellar dynamics 



1. Introduction 

The components of the Milky Way sample the diverse forms of matter in the universe 
from the black hole at the dynamical center to the dark matter component of the halo. 
Hypervelocity stars (HVSs), a class of recently discovered objects, are a probe linking the 
central black hole and its history to the outer halo. HVSs probably originate from interactions 
with the black hole; they travel out into the halo and beyond, serving as tracers of the 
Galactic potential on scales from 5 pc to 10 5 pc. 

To make testable predictions of the relationships between the observable properties of 
HVSs and aspects of the Galactic center and ha lo, we adop t the massive black hole (MBH)- 
binary star interaction first proposed by Hills (IHillsl Il988l ) as the injection mechanism for 
HVSs. We then explore the impact of the Galactic potential, the properties of binaries in 
the Galactic Center, the mass function in the Galactic Center, and stellar evolution on the 
observable spatial and velocity distributions of HVSs as a function of stellar mass. 



1.1. History and Challenges 



Hillsl (119881 ) showed that HVSs are inevitable when s ufficiently tight binary stars are 



disrupted by the tidal field of the central MBH. Much later, lYu fc Tremaind (120031 ) analyzed 
the production of HVSs and argued that the Hills production rates were excessive. They 
considered other mechanisms for the production of hypervelocity stars including interactions 
of single stars with a binary central black hole. 
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In 2005, iBrown et al.l (120051 ) reported the first discovery of a HVS. The star, SDSS 
J090745. 0+024507, is leaving the Galaxy with a heliocentric radial velocity of 853±12 km s _1 
and has the largest velocity ever observed in the Milky Way halo. Su bsequent photometr y 
revealed that the object is a slowly pulsating B main-sequence star (IFuentes et al.l 120061 1 . 
Only interaction with a massive black hole can plausibly accelerate a 3 M main-sequence 
B star to such an extreme velocity. 



The IBrown et al.l (120051) discovery ins pired further observations along with a wealth of 



Edelmann et all J2005h discovered a hypervelocity ma in se quence B star 



theoretical models. 

plausi bly ejected from the LMC (see also lBonanos &: Lopez-Moralesll2007l ) and lHirsch et al. 
(120051 ) reported a helium-rich subluminous O star probably ejected from Sgr A*. These three 
first discoveries were serendipitous. 



Brown et al.l (j2006al ) then carried out a targeted search for hypervelocity stars by using 



the SDSS to select faint B-type stars with lifetimes consistent with travel times from the 
Galactic center but that are not a normally expected stell ar halo population. S o far t his 
extensive survey has yielded another 7 hypervelocity stars ( Brown et al. 2006al lbl. l2007al lbl). 
The velocity of each of these stars exceeds the escape velocity from the Milky Way. 

The observational discoveries have inspired theorists to reconsider production mecha- 
nisms and ejection rates of hyperve locity stars . Ejection mechanisms include the origina l 



Hills binary encounter with a MBH (jHillslll988l ; lYu &: Tremaine 



encounters b e tween single s tars and binary black holes (lYu fc Tremaine 



2003 



Bromlev et al. 2006) 



2003 



Sesana et al. 



20061 . l2007al lbl:l Memttll2006l ). encount ers of single stars with an intermediate mass black hole 
(IMBH) in s piraling toward a MBH ([Hansen fc Milosavljeviq 120031 ; iGualandris et al.l 120051 ; 
Levinlbood : iBaumgardt et al.ll2006l ). and single star encount ers with a c luster of stellar mass 
black holes (SMBHs) around a MBH Jo'Learv & Loeblboosh . Recently Eu et all (bo07h pre- 
dicted hypervelocity ejection of tight binary stars in interaction with a central binary MBH. 
Discovery of a single tight binary HVS would be a strong indication of a central binary MBH. 

Absent the discovery of a binary HVS, distinguishing among the ejection mechanisms 
requires predictions of the resulting spatial and velocity distribution of the HVSs along with 



computatio n of expected ejection rates (e.g. Sesana et a l 2006, 2007a 



Hilla ( 119881 ) mechanism may produce the largest velocities dSesana et al. 



black hole may produce an anisotropic distribution of ejecta (jSesana et al. 



For example, the 



2007 af) a nd a binary 



2006|). 



Dynamical and evolutionary considerations for stellar populations near the Galactic 



Center also constrain HVS ejection mechanisms. For examp 
population of B stars at the Galactic Center favor the 



involving an inspiraling IMBH or a cluster of SMBHs (IPerets 



e, observations of the current 



Hills (1988) mechanism over models 



20070 . 
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Ejection rate estimates depend on an understanding of the way stars are scattered onto 
orbits intersecting the MBHs "loss cone" and on the specific ejection mechanism. Perets 
et al (2007a, 2007b) argue that scattering by giant molecular clouds gives HVS production 
rates from binary disruption by a single MBH that are consistent with the observations. 

In additio n to the HVSs, all o f the ejection mechan isms produce a bound population 
of ejecta (e.g., iBromley et al.ll2006l ). iBrown et al. j2007cM demonstrate the very probable 
presence of a bound population of ejected B stars. These stars have velocities between 275 
km s _1 and 450 km s _1 . In this velocity range the Brown et al surveys detect 18 outgoing 
stars and only 1 incoming star, consistent with the ~ 200-300 Myr lifet ime of these main 
sequence stars and a Galactic center origin (see also ISvensson et al.ll2008l ). 



Be cause HVSs travel across the Galaxy, th e y are powerful tracer s of the Galactic po- 
tential JGnedin et all 120051 : lYu fe Madaul l2007h . lYu fc Madaul (b()07t ) emphasize that the 
population of bound and unbound ejecta constrain the anisotropy of the halo independent 
of the ejection mechanism. 

The power of HVSs as probes of the Galactic center and as tracers of the Galactic 
potential provides strong motivation for acquisition of larger samples of these objects. Brown 
et al. focus on B stars because these objects are not an indigenous population in the halo 
and they are observable to large distances where the contrast between the density of HVSs 
and any indigenous population is largest. iKollmeier &: Gouldl (120071 ) suggest searches for 
possibly more numerous old-population, lower mass HVSs. Detection of lower mass HVSs 
could provide strong tests of models for their origin. 



1.2. Prom the Galactic Center to the Halo 



In the iHiha (119881 ) mechanism, HVSs naturally sample both the stellar mass function 
and the properties of the binary population in the Galactic Center. The ejecta then traverse 
the Galaxy and become observable in the halo. We show that the observable velocity and 
spatial distributions of HVSs are sensitive to the Galactic potential, particularly on scales 
of 5-200 pc from the Galactic Center. These distributions also reflect stellar evolutionary 
timescales. 

To include as complete a picture of the consequences of the Hills ejection mechanism 
as possible, we consider a broad range of issues in this paper. In §2.1 we simulate the 
MBH-binary intera c tion following; the procedures of IBromley et al.l (120061 ). Motivated by 



Kollmeier &: Gouldl (120071 ). we extend the simulations to cover lower mass stars and unequal 
mass binaries. In §2.2 we derive forms for the Galactic potential which are a good represen- 
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tation of the observations from 5-10 5 pc. Our results show that most of the deceleration of 
HVSs occurs within the central 200 pc of the Galaxy; thus, a proper fit to the observations 
in this region is a crucial step in understanding HVSs. 

In §2.3, §2.4, and §2.5, respectively, we describe our integration of HVS orbits through 
the Galaxy, the procedure we use to construct mock catalogs of HVSs, and the impact of 
azimuthal deflections on HVS orbits. 

We continue with more detailed discussions of the observable properties of HVSs in a 
simple spherically symmetric potential (§3) and in a three-component Galaxy (§4). One 
of the most striking general results of this investigation is that a central potential which 
matches observations in the central 200 pc of the Galaxy acts as a high pass filter preventing 
lower velocity HVSs from reaching the halo. 

In §3 and §4 we predict the spatial and velocity distributions of HVSs as a function of 
stellar mass and we address the underlying physical explanations for the behavior of these 
distributions. In §5 we connect the predictions with the observations by calculating relative 
space densities as a function of stellar mass (§5.1). This discussion emphasizes that the 
relative numbers of bound stars are sensitive to the mass function in the Galactic Center 
and relative stellar lifetimes; relative numbers of unbound stars are nearly independent of 
stellar lifetime. 

In §5.2 and §5.3 we analyze some sample HVSs search approaches. We argue that 
searches become increasingly difficult with decreasing stellar mass (§5.2). We also demon- 
strate that short lifetimes make the discovery of post main sequence HVS unlikely. We 
conclude in §6. 



2. The Simulations 



In iBromley et al.l (120061 ). we developed initial methods for predicting the observable 
velocity distribution of HVSs in the galactic halo. We derived ejection velocities appropriate 
for a single MBH disrupting binaries with 3-4 M primary stars and a realistic range of initial 
binary semimajor axes, au n . We then adopted a simple prescription for the Galactic mass 
distribution, integrated the orbits of stars ejected from the Galactic Center, and computed a 
first approximation to the full observable phase-space distribution of ejected stars on radial 
orbits in the halo. This analysis showed that the predicted median speeds of ejected stars as 
a function of distance in the halo are consistent with current observatio ns. We also p r edicted 
a population of ejected stars on bound radial orbits, then discovered by lBrown et al.l (j2007al ) 
as a 3.5a excess of B-type main sequence stars with velocities of 275-400 km s _1 . 
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Here, we extend our analysis to consider a wider range of binary star masses and sev- 
eral galactic potentials. Our goals are (i) to understand the relative frequencies of HVSs 
and bound ejected stars as a function of stellar mass and (ii) to explore the sensitivity of 
observable HVS radial velocities on the Galactic potential from 5 pc to 10 5 pc. 



As in iBromley et al.l (120061 ) , we simula te the ejection of stars from the Galactic Center 
in three steps. Following Hills! ( 1988 . 1992 ). we integrate orbits of binary stars passing by 
Sgr A* to quantify ejection probabilities and velocities. Then we use a Monte Carlo code 
based on semianalytical approximations for rapid generation of simulated catalogs of ejected 
stars in the Galactic Center. Finally we integrate the orbits of these stars in the Galactic 
potential to calculate how they populate the Galaxy's halo. 

To generate observable samples, the simulation incorporates the stellar main sequence 
lifetime from published stellar evolution calculations. We depart from lBromley et al.l (120061 ) 
by considering primary and secondary stars with masses mi and m2 in the range 0.6-4 M Q . 
In addition to considering orbits of stars in a simple spherical Galaxy, we also derive velocity 
distributions for stars ejected into a standard Galactic potential consisting of a bulge, disk, 
and extended halo. 



2.1. Massive Black Hole— Binary Star Interaction 



To obtain the spectrum of ejection speeds as a function of the proper ties of the dis 
rupted binary star, we first simulate the MBH-bi nary interaction. As in iBromlev et al. 



(2006 ), we use a sixth-ord er, symplectic integrator (jYoshidalll990l ; IBromley fc Kenyon 



2006 



Kenyon fc Bromleyll2006l ) to simulat e the disruption o f a binary system falling into an MBH 



with mass M. = 3.5 x 10 6 M Q (e.g. iGhez et al.ll2005l ). The binaries have a random orbital 
axis and phase; the semimajor axis is in the range aun = (flmm, a, max ), where a m i n depends 
on the binary mass and a max « 4 AU. We assume that the center of mass of the binary is 
on a hyperbolic orbit with an initial approach speed of 250 km s _1 (lHillslll988l ) and that the 
binary has negligible orbital eccentricity. 



Motivated by lKollmeier fc Gouldl (120071 ). we extend our simulations to equal mass bina- 
ries with lower mass primary stars and to unequal mass binaries (Table 1). Table 1 includes 
the main sequence lifetime t ms , the stellar radius at the time of core hydro gen exhaustion 



R ms , and the minimum au n for a binary with unit mass ratio, q = m 1 /m 2 = 1 (jSchaller et al. 



19921 ; ISchaerer et al.lll993l ). To minimize the probabili ty of a collision during the encounter 
10%), we set a min = 2 R ms ,i/rL (see also Table 1 of lGinsburg &: Loebll2007l ). where R ms ,\ 



(< 
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is the radius of the primary star and 



0.49g 2 / 3 



0.6g 2 / 3 + In (l+gV3) 



is the fractional radiu s of the inner Lagrangian surface for a binary system with a circular 
orbit (lEggletonlll983l ). This assumption allows us to sample a wide range of binary systems 
without worry that the interaction with the central bla ck hole leads to tidal in teraction 
between the binary components and possible coalescence (iGinsburg fc Loebl 120071 ). Table 1 
also lists a m i n for unequal mass binaries with 0.6 M Q secondary stars. 

Our calculations reprod uce published results and they are consistent with analytica l 
estimates of HVS production jHillJl988l . E992l : lYu & Tremaindl2~00~3l : bould & Quillenll2003h . 
The median ejection velocity - measured at infinity - is 



1760 



Qbin \ - 1 / 2 
0.1 M]) 



mi + rri2 
2 M 



1/3 



M, 



3.5 x 10 6 M m 



1/6 



f R km s" 



(2) 



where mi is the mass of the primary star and m<i is the mass of the secondary star. The 
factor fn is a normalization factor that depends on r c i ose , the distance of closest approach 
to the black hole: 



./ 



0.774 + (0.0204 + (-6.23 x 10~ 4 + (7.62 x 10~ 6 + 
(-4.24 x 10" 8 + 8.62 x 10~ n D)D)D)D)D, 



where 



D 



'"close 
^bin 



2M. 



1-1/3 



10 6 (m! + m 2 ) 



This factor also sets the probability for an ejection, P e f 

P ej w 1 - D/175 



(3) 



(4) 



(5) 



for < D < 175. For D > 175, r c i ose ^> a^; the binary does not get close enough to the 
black hole for an ejection and P e j = 0. 

For any binary system, the ejection probability and median ejection speed are set by 
the binary parameters - a^, m 1; and m 2 - and the distance of closest approach to the 
black hole r c \ ose . We assume that we can ignore the chance of a merger event between 
binary companions and that ejections of the primary and secondary star are equally likely. 
Because they can only occur in rare cases where a bin and r c i ose are both small (e.g., also 
Ginsburg fc Loebll2006l . 120071 ). mergers are unlikely to change our results by more than 10%. 
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Our simulations suggest that the binary mass ratio affects the ejection probability only 
when the mass ratio is extreme [mxl m i > 10). T here is then a consis tent preference for 
ejection of the secondary (for additional details, see lBromley et al.ll2006l ). When the binary 
contains unequal mass stars, the median ejection speeds of the primary and secondary stars 
are 



2m? 



mi + 1712 



1/2 



and 



V2 



2mi 



mi + rri2 



1/2 



(6) 



respectively. In all cases, the distribution of ejection speeds is approximately Gaussian with 
a dispersion of 20% of the median speed. 



2.2. The Galactic Potential 



Once ejected stars leave the central black hole, they travel out through the Galaxy 
and decelerate. Historically, mass models for the Galaxy developed for other applications 
are optimized to fit observations at distances r > 200 pc from the Galactic Center (e.g . , 



Dehnen fc Binnevlll998l : iKlvpin et al.ll2002l : lYu fc Tremainell2003l : IWidrow fc Dubinskill2005l ) . 



For ejected stars, however, the largest deceleration occurs in the central 200 pc (see Figure 
[1] below). Thus, we need a Galaxy model that reproduces the Galactic potential within 200 
pc of the Galactic Center. 

In thei r analysis of t he pro per motion of Arches, a cluster of young stars in the Galac- 
tic Center, IStolte et al.l (120071 ) also realized that standard Galactic potential models are 
inadequate to address dynamical phenomena in the Galactic Center. They derive a three 
component triaxial model for the central Galaxy which allows them to investigate the orbits 
of Arches cluster stars in the central potential. However, their application does not depend 
on an accurate potential beyond the Galactic Center; thus, their parameters do not match 
the potential at r > 0.5-1 kpc. 

Because HVSs traverse the Galaxy, we need a potential that fits observations from 
the Galactic Center to the outer halo. Thus, we derive Galaxy models that fit both the 
observations at 5-200 pc (where HVSs decelerate considerably) and at 10-100 kpc (where 
we observe HVSs today). 
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2.2.1. Parameterization of the Galactic potential 



We consider two parameterizations of the Galaxy potential. First, we follow iBromley et al 



( 120061 ) and adopt a simple, spherically symmetric density profile: 



p(r) 



C 



1 + (r/r c y 



To yield a mass of ~ 3 x 10 7 wi t hin 10 pc of the Galact i c Center (as inferred 
kinematics; Eckart fc Genzel 1997; Launhardt et al. 2002; Schodel et al. 2003; 



(7) 



rom stellar 



Ghez et al. 



2005) and a circular rotation speed of 220 km s -1 in the Solar neighborhood (IHogg et al. 
2005 ). we adoptQ r c = 8 pc and C = 1.396 x 10 4 M pc~ 3 . This density profile yields a mass 
of ~ 1.2 x 1 9 Mq inside 120 pc, c ompared to an estimate of 8 ± 2 x 10 8 M in the nuclear 
stellar disk (jLaunhardt et al. 2002 ). The corresponding potential is 



$ G (r) = -2nGCr 2 c [2 (r c /r) arctan(r/r c ) + ln(l + r 2 /r 2 c )}, (8) 

where G is the Gravitational constant. If we add the potential of the central black hole, the 
total potential is 

$( r ) = $ G (r) -GM./r (9) 

To examine the impact of $ on our results, we also consider a three component, spher- 
ically symmetri c potential that includes contributions from t he bulge ($&), disk (3>d), and 
halo (<&/i) (e.g.. iGnedin et al.ll2005l ; IWidrow fc Dubinskil 120051 1 . Specifically, we adopt 



G 



(10) 



where 



$ 6 (r) 

$d(R, z) 



-GM b /(r + a b ) 



-GM d /^Ri + ia d + (zi + b 2 y/^, 
-GM h \n(l + r/r s )/r 



and 



(e.g., iHernquistl Il990l ; iMiyamoto fc Nagail Il975l ; iNavarro et al.lll997l ). Here, r is the radial 
coordinate in a spherical coordinate system and (R, z) are cylindrical coordinates. 

To choose appropriate parameters for this potential, we consider multi-component mass 
models for the Galaxy. Most derivations adopt analytic functions for the mass density, 



1 Although Bromley et al. ( 20061 ) used the value for C quoted here in their numerical simulations, they 
quoted C = 1.27 x 10 4 in the text. Here, we quote the correct value for C. 
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p(R, z), and th en derive fitting parameters from least-squares fits to a set of dynamical ob 



serva b les (e.g., Dehnen &: Binnev 19981 ; iKlypin et all2002l ; iBattaglia et al.ll2005l ; iDehnen et al 



20061 ). IWidrow fc Dubinskil (120051 ) expand on this technique with an iterative approach that 
solves fo r a stable, self- consi stent bulge-disk-halo model which fits the standard observables 
(see also lWidrow et al.ll2008l ). Although both of these approaches fit the data well for Galac- 
tic Center distances ranging from ~ 200 pc to ~ 100-200 kpc , they are not optimized for 
the Galactic Center region at ~ 5-200 pc (see also lStolte et al.l 120071 ). 



Here, we derive a potential model that provides a reasonable m atch to obser v ations 

J2005h 



throughout the range 5-10 pc. Table 2 li sts our results . We follow iGnedin et al 



and adopt diskE and halo parameters from Klypin et al. ( 2002). For the 



parameters that fit velocity dispersion data at 5-200 pc (ITremaine et al. 



within 10 pc inferred from stellar kinematics flEckart fe Genzel 



Genzel et al. 2003: Schodel et al. 2003; Ghez et al 



disk inside 120 pc (8 ±2 x 10 8 M G 



Launhardt et al. 



200 



5) 



1997 



julge, we derive 
2002). the mass 



Launhardt et al. 2002 



and the mass of the nuclear stellar 



|2002| ). Using a Newton- Raphson method 
that over-weights measurements in the central 100 pc, we derive Mb = 3.55 x 10 9 M Q and 
rfe = 106 pc. 

Thi s model for the bul ge underestimates the stellar mass within a few pc of the Galactic 
Center ( IGenzel et al.ll2003l ). Because mass in the central 10-100 pc has the largest impact 
on HVS trajectories, this underestimate has little impact on our results. 

Table 2 compares the p arameters for our Gal axy model with other models us ed for calcu- 
lations of HVS trajectories. lYu fc Madaul (120071) ad opt the IKlypin et al.l (120021) parameters 
for th e bulge and disk and the iDiemand et al.l (120071 ) halo model. For the IDehnen fc Binnev 
( 119981 ) multi-component triaxial potential, we derived the best-fitting spherical, three com- 
ponent model for the radi al component of their pot ential. Our mass and scale length for the 
bulge agree well with the IDehnen fc Binneyl (119981 ) model. 



Table 3 compares the mass-radius relations for several Milky Way models with observa- 
tions. For the observed masses - listed as 'Target' in the Table - enc l osed w ithin 10 pc (M w ) 
and 120 pc (M 12 o), we adop t the kinematic mass fr om iGhez et al.l (120051 ) and the mass of 
the nuclear stellar disk from lLaunhardt et al.l (120021 ) . For the mass enclosed within r = 21 



pc (M 2 m), we use the measured velocity dispersion at 210 pc (136 km s 1 ; ITremaine et al. 



20021 ; IWidrow fc Dubinskil 120051 ) and make the conversion M210 = 3rv 2 / G. At larger radii , 



we adop t the circular r otatio n speed of 220 km s 1 at r = 8 kpc (for Msnnn; lHogg et al.ll2005l ) 



and the IKlypin et al.l (J2002|) Milky Way mass at 100 kpc (for M- 



100000 



Aside, from a low 



2 To maintain a circular velocity of 220 km s 1 at 8 kpc, we reduce the disk scale length from 5 kpc to 4 



kpc. 
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value at 210 pc, our mass model provides a good match to the target values. The lStolte et al. 



( 120071 ) results provide a b etter match at 210 pc, but underestimate the mass at 10 pc. Most 



other mass models (e.g., iDehnen fe Binney 



1998; Gnedin et a 



2005; 



Yu fe Madau 



2007 



Widrow fe Dubinskil (120051 ) 



) 



underestimate the mass substantially for all r < 210 pc. The 
models match the target data outside ~ 120 pc, but they overestimate the mass inside 100 
pc (see their Fig. 6). 



2.2.2. Impact of the central potential on HVS trajectories 



The radial variation of the acceleration, a = $/r, demonstrates the importance of the 
central 200 pc of the Galactic potential in the ejection of HVSs ( Figure ED). For r < 4 



pc, the central black hole dominates the deceleration of ejected stars (ITremaine et al.l 12002 



Widrow fe Dubinskill2005l ). At larger radii, the acceleration has a clear plateau with a radial 
width comparable to the scale length r b and then falls rapidly. For r > 5-10 r b , the disk 
and halo provide most of the deceleration. Although the disk adds to the halo component 
at r ~ 1-20 kpc, its contribution is always relatively small. 

To illustrate how the acceleration varies among different potentials, Figure [1] compares 
\a\ for our simple spherical potential with the three component model adopted here and 
a thr ee component mode l developed to analyze HVS trajectories in a triaxial Milky Way 



ilar to the 


Gnedin et al. ( 


2005 


) potential (e.g. 


Dehnen & ' 


3innev 


1998; ] 


Klvpin et al. 


2002 




Batta 


dia et al.ll2005: 


Dehnen et al. 


2006; 


Yu & Madaull2007 


) . Because the 


Widrow & Dubinski 



(120051 ) potential overestimates the mass in the Galactic Center, this model overestimates \a\ 
in the Galactic Center. 



The upper panel of Figure [T] compares \a\ for thelGnedin et al.l (120051 ) potential with our 
simple spherical model. The larger r b for the iGnedin et al.l (120051 ) potential yields a plateau 



with smaller acceleration (|a| ~ 5 x 10 7 cm s 2 ) at r ~ 20-200 pc. Other At r > 1-100 
kpc, the acceleration from the disk and halo match the simple potential. 

The lower panel of Figure [1] compares \a\ from our simple spherical model to \a\ for our 
adopted three component Galaxy. The relatively small scale length of our bulge yields a 
good match to the plateau (|a| ~ 6 x 10~ 6 cm s~ 2 ) of the spherical model at small radii, r ~ 
5-30 pc. The disk + halo potential matches the spherical model well at large distances, r > 
1-2 kpc. 



The acceleration profile of the inner Galaxy has a clear impact on the ejection of HVSs 
into the outer Galaxy (Figure [2]). The solid lines in Figure [2] show v(r), the variation 
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of velocity with Galactocentric distance for the simple spherical potential as a function of 
Vq, the initial ejection velocity at 1.45 pc. The fi gure also show s seve ral results for our 
three component Galaxy (dotted curves) and for the iGnedin et al.l (120051 ) potential (dashed 
curves). For the simple spherical potential, ejected stars with v$ < 625 km s -1 fail even 
to reach r ~ 1 kpc. Stars with 625 km s _1 < v o < 800 km s _1 reach r ~ 1-10 kpc; stars 
with vq > 800 km s _1 re ach r > 10 kpc. We derive similar results for our adopted three 
component Galaxy. In the lGnedin et al.l (120051 ) potential, stars with relatively small ejection 
velocities, 525 km s _1 < v o < 750 km s _1 can reach large halo distances, r ~ 10-100 kpc. 
Thus, the larger deceleration produced by a more realistic potential for the Galactic Center 
leads to smaller v (r) and less penetration of ejected stars into the halo. 

From these results, we exp ect clear diffe r ences between the observable velocity his- 
tograms of stars ejected into the IGnedin et al.l (120051 ) potential and the histograms of stars 
ejected into our three component Galaxy model. We discuss this issue in more detail in §3-4. 



2.3. Integration of Orbits through the Galaxy 



To generate populations of stars ejected from the Galactic Center, we perform Monte 
Carlo simulations of ~ 10 6 stars on radial orbits starting a small distance away from the 
black hole Sgr A*. Following iBromley et al.l ( 2006). we assume aun. is distributed with equal 



likeli hood per logarithmic interval (see also lAbtl 1 19831 ; iDuquennoy fc Mayorl 1 199 It iHeacox 
19981 ) between a m i n (Table 1) and a max ps 4 AU. Consistent with the expectation of strong 
gravitational focusing of the orbit of the center of mass of the binary by Sgr A*, we choose the 
distance of closest approach, r c i ose , from a distribution that varies as r close . For an adopted 
pair of stellar masses rti\ and m 2 , the random au n and r close yield the ejection probability, 
P e j. If P e j exceeds a random deviate, we choose an ejection velocity t> from a gaussian 
distribution with average velocity v e j and dispersion 0.2 v 



To follow the tr ajectories of these ej ected stars through the Galaxy, we use a simple 
leap-frog integrator (IBromley et al.l 120061 ) . We start each star on a radial trajectory with 
velocity v at a distance of 1.45 pc from Sgr A*. For our adopted potential, the mass in 
stars within this radius from Sgr A* is roughly 5% of the black hole's mass (Table [3]). Thus, 
our results are fairly independent of the starting position. Once the initial trajectory is set, 
we use the leap-frog integrator to derive the time-dependent velocity of the star through a 
Galaxy with potential $ (Eq. M or Eq. ITOl) . 



Simulations of HVSs with several potentials provide a better understanding of the phys- 
ical processes involved in ejecting stars from the Galactic Center into the Galactic halo. 
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Calculations with the simple spherical potential (Eq. [8]) allow us to isolate the impact of 
stellar evolution from the deceleration of ejected stars traveling through the Galaxy (§3). 
Simulations with the three component Galaxy (Eq. [TU]) yield clear estimates for the sensi- 
tivity of observable quantities on the mass in the Galactic Center and the triaxiality of the 
potential (§4). 



2.4. Predicting Observable Quantities 

To generate mock catalogs of HVSs, we assume the central black hole ejects stars at 
a constant rate on timescales comparable to the typical orbital timescales of the bound 
population, t or f, ~ 1 Gyr. We eject each star at a time t e j, randomly chosen between 
zero and the main sequence lifetime of the primary star t mSj i, and observe each star at a 
time t b s , randomly chosen between zero and the main sequence lifetime for the ejected star 
(t ms ,ej)- The first constraint (t e j < t ms> i) guarantees that the binary interacts with the central 
black hole before post-main sequence evolution disrupts the binary. The second constraint 
(t bs < t ms ,ej) guarantees that Earthbound observers detect the ejected star before it evolves 
off the main sequence. If t e j < t bsi we assign an ejection speed vq to the ejected star and 
integrate its orbit through the Galaxy for a travel time t t = t b s — t e j. For each ejected star 
with index i observed at time t bs,i ; the mock catalog includes the derived distance r i; the 
initial ejection velocity v^, the radial velocity at time t Q b s Vi, and the orbital elements for 
bound stars. 



Altho ugh this proced u re for generating catalogs of HVSs differs from the method de- 
scribed in iBromley et al.l (120061 ). it yields similar results. Moti vated by the high star 



formation rate i n young clusters near the Galactic Center (e.g., iFiger et al.l Il999l . 12002 



Lang et al.l 120051 ), we previously adopt ed a short timescale between formation and ejection 



e.g., t e j <C t mSt i (IBromley et al.l 120061 ) . Becau se binary stars ou tside nearby star forming 



regions can also encounter the black hole (e.g., iPerets et al.l 120071 ). here we adopt a less re- 



strictive constraint on the eject ion timescale (t e j < t ms ^). Direct comparisons of new catalogs 



with the IBromley et al.l (120061 ) catalogs show negligible differences in the predicted median 
speeds as a function of distance from the Galactic Center or in the velocity histograms of halo 
stars (see below). Thus, our simulations are relatively insensitive to assumptions concerning 
the ejection time from the Galactic Centejf] 



3 We assume a constan t Galactic potential. Because the Galaxy probably gains mass with time (e.g., 
Bullock fc Johnston! 120051 ). our estimates for the bound population of older, less massive stars are more 
uncertain than estimates for more massive stars. 
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2.5. Azimuthal Deflection of HVS Trajectories 



With main sequence lifetimes longer than a typical orbital timescale t or b ~ 1 Gyr, bound 
stars with masses smaller than 1.5 M can make several passes through the Galactic Center 
and, perhaps, interact with the central black hole or other stars near the Galactic Center 
before returning to the halo. For the spherically symmetric potential of Eq. we estimate 
that bound stars on trajectories that pass within 10 pc of Sgr A* have less than a 0.01% 
chance of passing within 10 AU of another star in the Galactic Center. The probability 
of encountering a star outside of the 10 pc region is even smaller. With speeds ~ 800- 
1000 km s^ 1 near the Galactic Center, a bound ejected star must pass within 0.1 AU of 
another star to suffer a significant deflection of its orbit. Thus , bound ejected sta rs do not 
change their orbits in a spherically symmetric Galaxy (see also lYu &: Madaul 120071 ). 



For the three component Galactic potential, significant deflections of HVS trajectories 
occur only if the potential is not spherical. Altho ugh a triaxial halo potential can produce 
significant proper motions in a sample of HVSs (IGnedin et al.ll2005l ). these produce little 
change in the radial compone nt of the velo c ity in the halo. After several passages through 
the Galactic Center, however, lYu fc Madaul (120071 ) note that a triaxial halo may redistribute 
a bound star's kinetic energy between the radial and angular velocity. For reasonable tri- 
axial halo potential s, the observed radial velocities of ejected bound stars change by ~ 10% 
JYu fc Madaul l2007h . 



A triaxial bulge potential can have a more dramatic impact on the observed velocities 
of ejected and bound HVSs . Current data suggest the mass distribution near the Galactic 



Center is not sph erical (see iDehnen &: Binneyi Il998l ; iKlypin et al.l 120021 ; iDehnen et al.l 12006 



Stolte et al.ll2007l . and references therein). Because it is easier to eject stars along the minor 
axis of a triaxial Galactic Center potential (e.g., Fig[TH3]), we expect potentially large dif- 
ferences in the observable velocities of HVSs as a function of galactic latitude. To estimate 
the magnitude of this effect, we derive predicted velocity distributions for ejected and bound 
stars in two different three component Galaxy potentials below (§4). 



3. Observable Ejected Stars in a Simple Spherically Symmetric Potential 

To isolate how the important physical effects of the ejection process relate to observable 
properties of ejected stars, we begin with a discussion of the evolution of ejected stars in 
a simple, spherically symmetric Galaxy (Eq. [7]). We describe the evolution in the three 
component Galactic potential in §4. We further restrict our analysis of observable properties 
to stars with r > 10 kpc, where current surveys can distinguish ejected stars from the 
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indigenous halo population and where the observe d radial velocity is ne ar ly indis tinguishable 
from the true radial velocity of ejected stars (e.g., Brown et al. 2006a b, 2007al lbl) . 



We divide these successful ejected stars into bound stars that orbit the Galactic Center 
and unbound stars that leave the Galaxy. Because our forms for the Galactic potential do 
not yield true escape velocities, we define unbound stars as ejected stars that reach r > 200 
kpc with v > +200 km s _1 (see Figure [2]). These stars have v > +450 km s -1 at r = 50 kpc 
and v > +300 km s _1 at r = 100 kpc. 

With typical orbital periods of 1 Gyr, bound stars falling back into the Galactic Center 
must have long lifetimes. Thus, bound stars are either low mass stars or massive stars on 
their first pass out through the G alactic halo. For un bound stars, the typical travel times 



to r ~ 40-80 kpc are ~ 100 Myr flBrown et al.l l2007bl ) . Thus, some 3-4 M Q unbound stars 



with lifetimes of 150-350 Myr evolve off the main sequence while traveling through the halo. 
This evolution removes potential HVSs from the catalogs, reducing the relative frequency of 
massive stars at r > 50-100 kpc (see below). Stellar evolution removes very few low mass 
stars from the catalogs; thus, the observable properties of unbound 0.6-2 M stars depend 
solely on the binary orbital parameters and the Galactic potential. 

Figure [3] shows the median ages for stars ejected from equal mass binaries on their 
first pass outward through the halo. Stars with radial speeds exceeding the escape velocity 
dominate the population of massive stars at all distances from the Galactic Center. Thus, the 
observed median ages of 3-4 M stars increase monotonically with r. Bound stars dominate 
the population of low mass stars at all Galactocentric distances. Thus, the median ages of 
0.6-1 M Q stars are roughly constant with r. For intermediate mass stars, the median ages 
represent a mix of mostly bound stars at small r and mostly unbound stars at large r. The 
distribution of median ages for 2 M Q stars is dominated by bound stars for r < r u « 20 kpc 
and by unbound stars for r > r u , where r u is the Galactocentric distance where the fraction 
of unbound stars exceeds 10%-15%. For 1.5 M stars, the bound population dominates for 
r ^ r u ~ 70-80 kpc. Thus, the median ages for 1.5-2 M stars fall for r < r u and rise or 
remain constant for r >r u . 

At large Galactocentric distances, the radial distribution of ejected stars depends on 
the stellar mass and the shape of the Galactic potential (Fig. |4J. Unbound stars (dashed 
lines in the Figure) are ejected at the largest possible velocities. They are the stars least 
affected by the Galactic potential and thus have shallow density profiles. For long-lived, 
low mass stars, the density profile of unbound stars is slightly steeper than the p oc r~ 2 
expected for a point-like potential. Bound stars are ejected at smaller velocities and reach 
smaller Galactocentric distances than unbound stars. Thus, bound stars have steeper density 
profiles. For both bound and unbound stars, low mass stars with main sequence lifetimes 
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longer than the travel time from the Galactic Center have shallow density profiles. Although 
they are typically ejected at larger velocities, many shorter-lived 3-4 M stars evolve off 
the main sequence before reaching 50-100 kpc and thus have the steepest density profile, 
particularly at large r. 

The large range in main sequence lifetimes among 0.6-4 M stars leads to clear dif- 
ferences in the fraction of unbound stars as a function of stellar mass and Galactocentric 
distance (Fig. Massive main sequence stars observed at r > 50 kpc must have relatively 
large ejection speeds and short travel times. Thus most massive stars at large r are unbound. 
Low mass main sequence stars ejected on bound trajectories have much longer to reach r > 
50 kpc before evolving off the main sequence. Because there are more ejected stars on bound 
orbits than on unbound orbits, most observable ejected low mass stars are bound at all r. 

The velocity distributions of observable ejected stars are also a strong function of stellar 
mass (Fig. [61 Table 2). For 3-4 M Q stars, the velocity histograms are asymmetric with 
median velocities of ~ 400-500 km s" 1 for stars with r > 30 kpc. The median velocity is 
sensitive to r: nearer stars have smaller median velocities. Because the typical travel time 
to r > 30 kpc is a significant fraction of t ms for a 3-4 M star, more massive stars must have 
larger ejection speeds to rea ch larger r. Thus, m ore distant regions of the halo tend to hold 



faster-moving ejected stars (IBromley et al.ll2006l ). 



For stars with masses < 2 M Q , the velocity histograms are more symmetric with de- 
creasing stellar mass. Two physical effects produce symmetric velocity histograms for the 
lowest mass ejected stars: (i) most ejected stars that reach r > 10 kpc are bound and (ii) 
because 1 M stars live far longer than 2-4 M stars, most of the bound population consists 
of low mass stars (Fig. [5]). For typical t or b ~ 1 Gyr, these stars live long enough to make 
several orbits of the Galaxy. Bound 2-4 M stars nearly always evolve off the main sequence 
before completing a single orbit of the Galaxy. Thus, bound stars increasingly dominate the 
velocity histograms of observable lower mass stars. 

To illustrate the importance of the stellar lifetime in more detail, Figs. [7H8] show the 
median speed as a function of Galactocentric distance for the entire observable population 
(Fig. [7]) and for outgoing observable stars only (Fig. [S]). The 3-4 M stars ejected at 
relatively small velocities (vq < 700-800 km s _1 ; Fig. [2]) do not live long enough to reach 
r > 10 kpc. Lower mass stars live long enough for nearly all ejected stars to reach r > 10 
kpc. Thus, more massive stars observed at r > 10 kpc move faster than lower mass stars. 
At even larger distances, the short stellar lifetimes of massive stars select for stars ejected at 
even larger velocities. Thus, the median sp eed of observable massive ejected stars is a strong 



function of r (see also lBromley et al.l 120061 ) . The longer lifetimes of low mass stars allows a 



larger fraction of ejected stars to reach large r; thus, the median speed depends only weakly 
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on r. 

The median speeds of primary stars ejected from unequal mass binaries are fairly in- 
sensitive to the mass ratio q of the disrupted binary (Fig. [§]). Massive primary stars ejected 
from low mass binaries tend to have somewhat smaller median speeds than massive stars 
ejected from more massive binaries. For binaries with 0.6 M secondaries, the large mass 
ratio leads to the smallest ejection velocities and thus the longest travel times for the 4 M 
primary (Eq. ED- For binaries with 1 M secondaries, the longer travel times lead to smaller 
median speeds only for r > 80 kpc, where the density profile is very steep (Fig. H]). 

The median speeds of secondary stars are much more sensitive to q (Fig. [TO]) . Although 
secondary stars ejected from unequal mass binaries are ejected at larger speeds than their 
equal mass counterparts (Eq. EJ), two features of the ejection process lead to smaller median 
speeds for secondary stars observed at r ~ 10-80 kpc (compare with Fig. [7]). Because 
the primary stars of unequal mass binaries have shorter lifetimes than the secondaries, the 
median age of an ejected secondary star is smaller for an unequal mass binary than for 
an equal mass binary. Thus, the secondary stars in unequal mass binaries have longer 
median travel times. Longer travel times allow secondaries ejected at velocities larger than 
the median ejection velocity (Eq. [2]) to reach r > 80 kpc, removing these stars from the 
population at r ~ 20-80 kpc. Longer travel times also allow more slowly moving ejected 
stars to reach large r. Removing the highest velocity stars from the population and adding 
more slowly moving stars to the population reduces the median velocities of secondaries at 
r ~ 20-80 kpc. 

To summarize, simulations in the simple spherical potential yield two classes of ejected 
stars observable at r > 20 kpc. The bound population of mostly long-lived, low mass stars 
has a steep radial density profile, p(r) oc r~ n with n m 3, and a symmetric radial velocity 
profile centered at v r ps km s _1 (Figs. 0H2J Table BJ. Short-lived 3-4 M Q bound stars have 
slightly steeper radial density profiles and provide a high-velocity tail to the radial velocity 
distribution at +275 to +500 km s~ l . Unbound ejected stars with v > +500 km s" 1 have 
shallow radial density profiles (n ~ 2-2.5). The short lifetimes of 3-4 M Q unbound stars 
yield very steep density profiles at large radii {n > 3 for r > 80 kpc). Thus, the relative 
frequency of 3-4 M unbound stars decreases relative to lower mass unbound stars at r > 
80 kpc. 
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4. Observable Ejected Stars in a Three-Component Galaxy 

We now consider the evolution of observable ejected stars in a three compo nent potential 



consis ting of a bulge, disk, and halo (Eq. [10]). We discuss results for the iGnedin et al. 



(120051 ) Galaxy model, which has a relatively small acceleration at small r, and for our three 
component model Galaxy, which has a larger acceleration at small r. The acceleration 
at small r for our model Galaxy is comparable to the acceleration for the simple spherical 
potential (§2.2; Fig. [I]). To minimize differences in median velocities caused by the relatively 
massive disk in our model Galaxy, we derive results for ejections a long the z ax i s. To 



minimize the number of low velocity stars ejected into the halo of the IGnedin et al.l (120051 ) 
model Galaxy, we derive results for ejections in the disk plane, where the larger deceleration 
from the disk prevents more low velocity ejected stars from reaching 10 kpc. 

At 10-60 kpc, the velocity histograms derived for ejections into our model Galaxy and 
the simple spherical Galaxy are indistinguishable (Figure fTTi) . In both cases, the distributions 
for 4 M Q stars are greatly skewed to large velocities, with more low velocity stars at 10-30 
kpc than at 30-60 kpc. These models yield a small (~ 10%) fraction of stars returning to 
the Galactic Center with negative velocities. For 0.6 M stars, the distributions are nearly 
symmetric about zero velocity, consistent with a large population of bound stars in both 
examples. In addition to a modest fraction (~ 10%) of unbound stars, models for 0.6 M Q 
stars have a large fraction (~ 45%) of returning stars. 



At 10-60 kpc, calculations with the IGnedin et al.l (120051 ) Galaxy potential yield more 



symmetric velocity histograms and many fewer unbound stars than calculations using the 
other Galaxy models. For 2-4 M stars, the velocity histograms are skewed to high velocities, 
but have a smaller fraction of stars with velocities exceeding 600-800 km s -1 (Figure [TT]) . 
Ejections into this potential also produce a much larger population of stars at 10-30 kpc 
than at 30-60 kpc. For lower mass stars, the velocity histograms are more similar to those 
derived from the other Galaxy models. However, velocity histograms with this potential 
have smaller dispersions about the median velocity. 

Table H] compares median speeds at 30-60 kpc for the three model Galaxies. For 2-4 M Q 
primary and secondary stars, the median speeds derived for our three component Galaxy 
are ~ 5% to 15% smaller than those derived for the simple spherical Galaxy. This difference 
results from the potential at r ~ 1 kpc, where the three component potential produces a 
smalle r acceleration t han t he simple potential (Fig. [T]). The median speeds of massive stars 



in the IGnedin et al.1 (120051 ) potential are typically ~ 75% of those for the other potentials. 



For 0.6-1 M stars, the median speeds derive d for the thre e mode ls are similar. Although 



the median velocities of low mass stars in the IGnedin et al.l (120051 ) potential are ~ 2/3 the 



median velocities derived for the other potentials, the typical difference of 6-12 km s 1 is 
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small compared to current observational errors. 



In these calculations, the large differences between the median speeds and the shapes of 
the velocity histograms are set by the Milky Way potential at r = 5-200 pc (Figured]). At 
small r, the potential acts as a high pass f il ter, p reventing low velocity stars from reaching 
the halo (Figure [2]). In the iGnedin et all (120051 ) potential, lower velocity stars reach the 
halo. Thus, the median speeds of observable stars are smaller and the velocity histograms 
are narrower. In the simple spherical Galaxy and our three component Galaxy, low velocity 
stars are trapped at r < 10 kpc. Thus, the median speeds are larger at all r > 10-20 kpc. 
Because low velocity bound stars remain at r < 10 kpc in our three component Galaxy, the 
velocity histograms of these stars are also broader for all r > 10-20 kpc. 

Because the potential of the central 200 pc of the Galaxy makes such a large difference 
in the median speeds of observable HVSs at 10-60 kpc, our results suggest that HVSs might 
be a useful probe of the triaxial potential of the cen tral Galaxy. To te st this possibility, 
we consider tw o triaxial models f or the inner Galaxy. IStolte et al.1 (120071 ) analyze the mass 
density data of lLaunhardt et al.1 (120021 ) and derive a three component triaxial model for the 
central 200 pc. Each component in this model has axial ratios Xq : jjq : z Q . With a minimum 
axial ratio of 0.71 in the z-direction (for the nuclear stellar disk), we predict an expected 
reduction of 15-20 km s _1 in the median speeds of stars ejected along th e z-axis relative to 



stars e jected in the plane of the Galaxy. Using OGLE data for the bulge, iRattenbury et al. 
( 120071 ) derive a triaxial model for the Galactic bar with x : Uo '■ z o = 1-0 : 0.35 : 0.26 



Assuming that the rotation of the bar eliminates any signal in HVS speeds along the y-axis, 
we expect a reduction of 25-30 km s -1 in the median speeds of HVSs ejected along the z-axis 
of the bar relative to those ejected in the Galactic plane. 

Given current uncertainties in mass models for the Galaxy, our simulations suggest that 
the radial density profiles and velocity histograms of HVSs at r > 20 kpc are most sensitive 
to the form of the Galactic potential for r < 0.1-1 kpc. At r > 1 kpc, small changes 
to the potential produce negligible differences in observable quantities. At r < 200 pc, 
triaxial potentials consistent with observations produce ± 10-30 km s -1 differences in the 
median velocities of observable unbound stars as a function of galactic latitude and negligible 
differences in the radial density profiles. These differences are comparable to the variation 
in median speeds among 0.6-1.5 M ejected stars, where the shorter stellar lifetimes of more 
massive stars produce median speeds offset from v r ~ 0. With high quality stellar mass 
estimates, HVSs could provide constraints on the triaxiality of the central galactic potential. 
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5. Scaled Density Profiles and Searches for HVS 



The K ills ( 



9881) ejection mechanism for HVSs has several observable consequences. In 



Bromley et al.l (120061 ) and in this paper, we show that radial velocity surveys for massive 



main sequence stars should reveal comparable numbers of bound and unbound HVSs. Our 
results also suggest that 4 M HVSs should have smaller space densities at r ~ 50-100 kpc 
than 3 M HVS. Observations confirm both predictions. Brown et al. ( 2007al fbh identify 26 
bound B-type HVSs with v r > 275 km s -1 and only 1 bound HVS with v r < —275 km s -1 , 
supporting the notion that the B-type HVSs are main sequence stars and are not horizontal 
branch stars. In these data, longer-lived 3 M Q HVSs fill the survey volume out to r ~ 
80-90 kpc; shorter-lived 4 M HVSs are missing at large distances. This result confirms our 
expectation that stellar lifetimes set the HVS space density at large r. 



From an analytic analysis of the radial density profiles of unbound stars, iKollmeier fc Gould 



( 120071 ) develop a strong motivation to search for HVSs among 0.6-1 M Q stars. They con- 
clude that a survey of faint stars near the main sequence turnoff might yield as many as 
0.022 HVSs deg -2 , roughly a factor of 20 larger than the ~ 0.001 HVSs deg -2 inferred from 
our targeted survey of B-type stars (IBrown et al.ll2007bl ). Our calculations in §3-4 suggest 
that the large bound population of observable ejected turnoff stars overwhelms the unbound 
component, especially at r ~ 10-40 kpc where the radial velocities of turnoff stars can be 
measured with existing multi-objects spectrographs (Fig. E]). However, if the relative space 
density of bound, low mass ejected stars is large enough relative to the i ndigenous halo pop 



Kollmeier fe Gould 



Kollmeier fc Gould 



ulatio n, fruitful searches for low mass HVSs might still be possible (see 
20071) . The success of the HVS ejection model for 3-4 M stars and the 
(120071 ) suggestion leads us to take a more detailed look at the likelihood of identifying large 
samples of HVSs among lower mass stars. 



5.1. Relative Space Densities as a Function of Stellar Mass 

To quantify the relative space densities of HVSs as a function of stellar mass in our 
simulations, we must scale the density profiles in Fig. |H For each stellar mass, we require 
the relative frequency of primary and secondary stars as a function of stellar mass (the mass 
function), the relative frequency of binaries with the range of separations (a m i n , a max ), and 
an efficiency factor relating the number of stars ejected from the Galactic Center to the 
number observed at r « 10-100 kpc. For simplicity, we assume that the binary frequency is 
independent of primary mass and that the mass function is independent of orbital separ ation. 



Altho ugh observations suggest the binary frequency depends on stellar mass (e.g., lLada 



20061 ). the differences in binary frequencies among nearby 0.6-4 M stars are small compared 
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to uncertainties in the mass function and the relative frequency of binary separations. 



To construct a first estimate of the scaled density profiles of unbound stars as a function 
of stellar mass, we develop a simple model. For each stellar mass m, we derive a weighting 
factor 

x(m) = xx(m) x 2 (m,a) x 3 (m) t ms (12) 

where x\{m) is the relative number of stars of mass m, X2(m,a) is the relative number 
of binary stars with orbital separations in the range a^n = (amin^max), and x^{m) is the 
relative fraction of stars ejected from 1.45 pc that reach r > 10 kpc. To account for stellar 
evolution, we scale the weighting factor by the main sequence lifetime t ms in Table [TJ We 
set x(m = 0.6 M ) = 1 and normalize other weights accordingly. This model assumes that 
primary and secondary stars are selected from the same mass function, a reasonable first 
approximation. 



To derive X\, we adopt a simple power law mass function 

Salpeter 19551 ; Miller fc Scalo 1979 ). Integrating this functi 



(13) 



with q « 1-1.5 (e.g. ISalpeterl 119551 : iMiller fc Scald 119791 ). Integrating this function over m 
yields relative numbers of stars in mass ranges with lower mass limits mi and upper mass 
limits m u . For a set of stellar masses, S{m) = {0.6,1.0,1.5,2.0,3.0,4.0}, we choose a set 
of lower mass limits S(mi) = {0.4,0.8,1.25,1.75,2.5,3.5} and a set of upper mass limits 
S(m u ) = {0.8, 1.25, 1.75, 2.5, 3.5, 4.5}. 

To guide our choices for the exponent q, we re l y on recent obse rvations of the Arches 
cluster near the Galactic Center (IStolte et al.ll2005l ; iKim et al.l 120061 ). These analyses sug- 
gest a current mass function with q « 0.9, top-heavy compared to the local Salpeter mass 
function with n m 1.35. However, dynamical interactions among cluster stars may produce 
mass segregation on timescales comparable to the cluster age, flattening th e mass function 
in the core and steepening the mass function in the outskirts of the cluster (IDib et al.l 120071 : 
Portegies Zwart et al.ll2007l ). Thus, the most likely exponent for Arches stars is probably q « 
1.0-1.35. To provide estimates for a range of q, we derive S(xi) = {1., 0.29, 0.13, 0.086, 0.05, 0.025} 
for q = 1.35 and S(x{) = {1., 0.36, 0.18, 0.14, 0.09, 0.05} for q = 1.0. 

Observations of the Arches cluster and other young stars currently provide no con- 
straints on the frequency of separations for close binaries in the Galactic Center. To de- 
rive X2, we thus rely on observations in the local solar neighborhood. For nearby binaries 
with A-type and B-type primary s t ars, au v is d istributed with roughly equal likelihood 
per logarithmic interval (lAbtl Il983l : iHeacoxl Il998l ). For stars with orbital separations in 
the range au n = (amin,a m ax) from Table 1, this assumption yields a set of values for x%, 
S^xz) = {1., 0.66, 0.45, 0.41, 0.39, 0.39}. 
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For nearby sola r-type stars, the frequency of binary separations is closer to a log- 



normal distribution (jDuquennoy &: Mayor! Il99ll ). For this relation, we derive S^a^ 



{1., 0.95, 0.90, 0.87, 0.86, 0.85}. Because this frequency distribution under-weights binaries 
with small separations, it under-weights low mass ejected stars relative to massive ejected 
stars. To maximize the predicted relative density of low mass stars, we adopt the Si set of 
weights for X2- This decision may overestimate the relative abundance of solar- type starg^l 
by a factor of ~ 2 

To derive 23, we rely on the calculations for §3-4. For a complete derivation of x%, we 
require a full set of simulations spanning the entire range of binary separations and mass 
functions for primary and secondary stars. Because very few binaries produce ejected stars, 
this simulation requires a significant computational effort for little return. In the set of sim- 
ulations made for this paper, the ratio of failed to successful ejections is nearly independent 
of the masses of the primary and secondary stars and the initial binary separation. Thus, 
for our first estimate of x, we adopt x 3 = 1. 

Fig. [T21 shows the scaled relative density profiles of observable main sequence stars from 
this simple model. For clarity, we normalize the density profiles for 1-4 M stars to the 
density profiles for 0.6 M stars. Thus, the relative densities of 0.6 M stars are unity at 
all r. For all 0.6-4 M ejected stars (upper panels), the density profiles scale inversely with 
stellar mass. Thus, low mass ejected stars are relatively more abundant than massive stars at 
all r. Because shallower mass functions contain relatively more massive stars, the differences 
in the scaled density profiles are smaller for smaller n. For unbound 1-4 M stars at r ~ 
10-200 kpc, the density profiles are roughly constant with r and scale inversely with stellar 
mass (lower profiles). Aside from the dramatic decline in the relative density of 4 M stars 
at r > 30 kpc due to their relatively short main sequence lifetimes, massive unbound stars 
are a factor of ~ 10 less abundant than low mass unbound stars at all r. 

Two physical effects cause the changes in the slopes of the scaled density profiles in Fig. 
H2l Because massive ejected stars contain a larger fraction of unbound stars (Fig. Hj), massive 
ejected stars have shallower density profiles at small r than low mass ejected stars. Thus, 
the scaled density profiles for massive stars rise with r. At large r, however, short stellar 
lifetimes remove massive stars from the sample of observable unbound stars. For 4 M stars, 
the typical travel time to r ~ 50-100 kpc is comparable to the main sequence lifetime, t ms ~ 
160 Myr (Table 1). Thus, the relative densities of 4 M stars decline considerably at r ~ 



4 Other physical processes - includin g binary evaporati on near the MBH (jPereta 120071 ) and interactions 
between binaries and molecular clouds (jPerets et al.M2007T) - can lower the relative abundance of solar-type 
stars among HVS. To maximize the predicted density of low mass stars among HVS, we ignore these processes 
in our estimate of x-i. 
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50-100 kpc. For 2-3 M stars, the typical travel time to r ~ 50-100 kpc is less than half of 
the main sequence lifetime, t ms « 300 Myr to 1 Gyr. Thus, the relative densities for these 
stars decline with r less dramatically than the relative densities for massive stars. 

The relative numbers of bound stars as a function of stellar mass depend primarily on 
the mass function and the relative stellar lifetimes (Eq. (fl2|) ). For a mass function with q = 
1.35 and the stellar lifetimes in Table (H the relative number of 0.6 M and 3 M Q stars is ~ 
800:1, close to the prediction of ~ 1000:1 from the detailed model. 

In contrast, the relative numbers of unbound stars are nearly independent of stellar 
lifetime. Because the time for unbound stars to reach r ~ 50-100 kpc is a small fraction of 
the main sequence lifetimes for 0.6-2 M Q stars, the relative numbers of 0.6-2 M Q unbound 
stars depends primarily on the adopted mass function. Thus for q = 1.35, the relative number 
of 0.6 M and 2 M Q unbound stars is ~ 12:1. For more massive stars, the travel time is 
comparable to the main sequence lifetime. If we assume that ~ 50% of 3-4 M stars evolve 
off the main sequence as they travel from 10 kpc to 100 kpc, the relative numbers of 0.6 
M and 3-4 M unbound stars is ~ 50:1 for a mass function with q = 1.35, close to the 
prediction in Fig. 



5.2. Predicted Detections 



Converting the relative density pr ofiles of Fig. [121 into predicted detection rates re- 
quires a targeting strategy. We follow iBrown et al.l (j2006al ) and assume target stars are 
selected from a magnitude-limited survey. For simplicity, we consider a shallow magnitude- 
limited survey to V max = 21 and a deeper survey to V max = 23. To estimate relative 
detection rate s , we derive absolute V magnitu des as a function of stellar mass using the 
Schaller et al.l (119921) and ISchaerer et al.l (119931 ) evolutionary tracks and bolometric correc- 



tions from iKenyon fc Hartmannl (119951 ). We use these estimates to derive the maximum 
observable distance, r max (m, V max ), for each stellar mass in our set of masses. Integrating 
the relative densities in Fig. [T2]over distance from r = 8 kpc to r max (m, V max ) then yields 
the relative numbers of stars of each mass as a function of V max . For surveys that yield 
one observable 4 M HVS, we predict n e)Tn the relative number of observable ejected stars 
at mass m and n Ujm the number of observable unbound stars at mass m. Table [5] lists our 
results. 

Shallow surveys with V max < 21 are sensitive to HVSs of all masses (see columns (3) and 
(5) of Table [5]). Roughly 40% (12%) of the unbound (bound) ejected stars are 2-4 M stars 
with My = +1.7 to —1. Because lower mass stars have longer main sequence lifetimes, these 
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stars comprise a larger fraction of bound stars in a magnitude-limited survey. Nearly all (~ 
94%) observable ejected stars with m « 0.6-1.5 M are bound to the Galaxy. However, ~ 
25% of ejected 2-4 M stars are unbound. Thus, shallow surveys targeting 2-4 M stars 
have a higher probability of observing unbound HVS than surveys for lower mass stars. 

Deeper surveys with V max = 23 yield larger fractions of lower mass stars in the bound 
and unbound populations (see columns (4) and (6) of Table [5]). Although our analysis 
predicts that ~ 50% (75%) of the unbound (bound) ejected stars with V < 23 are 0.6-1 M 
stars, only ~ 5% of observable low mass ejected stars are unbound. Among 1.5-4 M stars, 
we derive a much larger fraction of unbound stars, ~ 25%. Because 3-4 M unbound stars 
are relatively rare in surveys with V max = 23, deep surveys targeting 1.5-2 M stars have 
the highest probability of detecting unbound HVS. 

Several competing physical effects combine to produce the results in Table [5J Mas- 
sive stars are much brighter than low mass stars; thus, magnitude-limited surveys probe 
much larger volumes for massive stars than for low mass stars. For typical density profiles 
p oc r~ 2 , the relative nu mbers of stars scale with their limiting distances r max (see also 



Kollmeier fc Gouldl 120071 ). All unbound stars spend comparable times traveling from 10 kpc 
to 200 kpc. Thus, the relative numbers of unbound stars scale with the mass function at 
the Galactic Center, r max , and the ratio of travel time to stellar lifetime. For 1 M stars 
and 3 M stars, these factors yield relative numbers of 6:1 (mass function), 1:10 (r max ), and 
2:1 (stellar lifetimes); the simple estimate for the relative number of ~ 1:1 is close to the 
tabulated values. 

To make a simple estimate for the relative numbers of bound stars, we need a larger 
correction for relative stellar lifetimes. Because low mass stars live longer than the typical 
orbital time of ~ 1 Gyr, the bound population samples low mass stars ejected throughout 
the lifetime of the Galaxy. For 1 M and 3 M stars, the relative stellar lifetimes (~ 30:1) 
yield relative numbers of 15:1 in a magnitude lim ited survey, close to the tabulated values 



of 10-20:1. Thus, as iKollmeier fc Gouldl (120071 ) first pointed out, low mass bound stars 



dominate the population of observable ejected stars in all magnitude limited surveys. 

As a final piece of this analysis, we conside r identifying HVSs within the large population 



of indigenous halo stars. (IBrown et al.ll2006al ) use SDSS (g 1 — r') and (u 1 — g')o colors to 



target 3-4 M HVS with g' Q < 19.5. Searching for massive stars in this color-magnitude 
space has several advantages over searches for lower mass HVSs: (1) surveys for bright, 
massive HVSs sample a large volume, (2) the predicted fraction of unbound stars is ~ 50% 
(Fig. H]), (3) the predicted velocity histogram for bound stars is more asymmetric (Fig. [6]) 
and thus they are more readily detectable, and (4) two colors reduce white dwarf and quasar 
contamination considerably, yielding a 3 M HVS candidate list of manageable size (see also 
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Brown et alJl2007al lbf) 



In a 6800 square degree survey of 575 candidates with 17.0 < g' < 19.5, iBrown et al. 



(j2007bl ) identify 7 unbound HVSs and 9 probable bound HVSs. The surface number density 
of unbound 3 M HVS is ~ 0.001 per square degree; the surface number density of observed 
3 M ejected stars is ~ 0.002 per square degree. In this shallow s urvey, the relative nu mbers 
of bound and unbound HVS agree well with model predictions (IBrown et al.ll2007bl ). The 
survey detection efficiency (number of HVS detected / number of candidates) is 1.2% for 
unbound objects. 



At the other end of the mass range we explore in Table [SJ iKollmeier &: Gouldl (120071 ) 
propose observing fainter (19.5 < g' < 21.5) stars near the main sequence t urnoff with 
masses m ~ 0.8-1.3 M^. A pplying scaling arguments to the lBrown et al.l (j2006al ) detections, 



Kollmeier &: Gouldl (120071 ) predict 1 turnoff HVS per 45 square degrees within their suggeste d 
survey limits. These arguments, and a later corroborating estimate by (IBrown et al.ll2007bl ). 
assume that the velocity distribution of HVSs is independent of their mass, contrary to the 
more detailed calculations for Table 



However, IKollmeier fc Gouldl (120071 ) scale the relative numbers of massive HVSs by the 
present day mass function, which is a product of the initial mass function (IMF) and stellar 
evolution^. This approach is valid for t he bound population. Thus, the predictions of 1 
turnoff HVS per 45-50 square degrees (IKollmeier fc Gouldl 120071 : IBrown et al.l l2007bl ) are 
reasonable estimates for the bound population of low mass HVSs. Because all unbound 
stars spend roughly equal amounts of time traveling through the halo, the relative numbers 
of unbound stars scale with the IMF. This scaling reduces the predicted sky surface density 
of unbound low mass HVSs relative to unbound massive HVSs as outlined above. Thus, 



Kollmeier fc Gouldl (120071 ) overestimate the detectable population of unbound HVSs by a 
factor of ~ 30-100. Unbound stars are a small fraction of the observable ejected turnoff 
stars. 

According to Table ~ 5% of turnoff stars in a survey for low mass HVS with g' < 
21.5 are unbound. Within the more dominant bound population, ~ 20% should have |i>| > 
275 km s -1 . For a predicted sky surface density of 1 bound HVS per 45-50 square degrees, 
detecting unbound or bound low mass HVS within a larger (300 per square degree) popula- 
tion of indigenous halo stars may be challenging. Detecting a single unbound low mass HVS 
requires observations of ~ 3 x 10 5 stars. This large sample would yield ~ 5 bound candi- 
dates with \v | > 275 km s™ 1 . We emphasize that only the signature of an unbound HVS is 
unambiguous; the bound component may be difficult to distinguish from the well-populated 



J Formally, we derive a crude estimate of the present day relative mass function as x\{m) t ms in Eq. ((T 
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tails of the velocity distribution of the indigenous population. 



These estimates are sensitive to the ejection model for HVS and to the adopted mass 
function and the fraction of stars in binaries with separations (a m i n , a max ) in our model (x\ 
and x-i in Eq. (fl2|) ). As outlined above, both of these inputs to our model are uncertain. 



Nonetheless, these models indicate that identification of HVS is increasingly difficult as 
the mass of the star decreases. The drivers of this conclusion are that (1) in a magnitude 
limited survey massive stars are visible to greater depth, (2) the fraction of unbound stars 
increases with mass, (3) the velocity distributions are increasingly symmetric for less massive 
HVS, and (4) the indigenous contamination increases with decreasing stellar mass. 

If HVS are ejected by a bina ry MBH instead of a sing le MBH, we exp ect d i fferent 



scaling laws with stellar mass (e.g., ISesana et al.ll2006l . l2007al ). Clearly, for the IHillsi (119881 ) 



mechanism, the properties of binaries in the Galactic center are important; they are irrelevant 



for ejection by a binary MBH. 



the origin of HVS (e.g., iPerets 



hus, searches for low mass HVS provide tests of basic inputs 



20071 ). 



to our estimates for the jHillsl (1988) ejections mechanism as well as the overall picture for 



5.3. Post Main Sequence HVS 



Although our results indicate that low mass main-sequence stars are not prime targets 
for dedicated HVS surveys, we must consider whether their post-main sequence descendants 
- horizontal branch (HB) and red giant branch (RGB) stars - are more favorable targets. 
HB stars are as luminous as 2-4 M main se quence stars and ar e thus observable to larger 
distances than their 1-2 M progenitors (see lBrown et al.ll2007bl ). With larger luminosities 
than HB stars, RGB stars are also reasonable targets for HVS surveys. 

To quantify the probabilities for detecting HB or RGB HVSs, we derive r ma Jm, V) from 



adop ted absolute magnitudes and lifetimes. For HB star s, we adopt My « 0.5 (IBrown et al. 



20081 ) and lifetimes ~ 1% of the main sequence lifetime (jYi et al.ll200ll ). Scaling the relative 
numbers by the relative lifetimes (a factor of 0.01) and the relative volumes (a factor of 7-8), 
the HB descendants of 1-2 M main sequence stars are a factor of ~ 10 less abundant than 
their main sequence progenitors. Thus, HB stars are unlikely to provide many HVSs in a 
deep survey. 

For RGB stars, we consider deep optical and near-infrared (IR) sur veys. We adopt 
an ab solute brightness at the tip of th e RGB (My £a —1.5; ~ —6; iBellazzini et al. 
20011 ) and a typical lifetime of 10 Myr (jYi et al.ll200ll ). For optical surveys with limiting 
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magnitudes V < 21-23, we predict that RGB stars are ~ 2% as abu ndant as 1-2 Mp mai n 
sequence stars. In an IR survey with a depth of K = 15 (e.g., 2MASS lSkrutskie et al.ll2006l ). 
we predict that RGB stars are < 1% as abundant as 1-2 M Q main sequence stars in a deep 
optical survey. Although RGB stars are observable to distances of 100-200 kpc in deep IR 
and optical surveys, their short lifetimes preclude detection as HVSs. 



6. Conclusions 

HVSs are a fascinating newly discovered class of objects because they connect the Galac- 
tic C e nter w ith the outer halo of the Milky Way. We explore these connections by using the 



Hilla ( 119881 ) model to inject stars into the Galactic potential. We track the journeys of these 



ejected stars across the Galaxy and derive simulated catalogs of observable HVSs. 

The foundation for our model includes the construction of forms for the Galactic poten- 
tial which fit observations over the range 5-10 5 pc. We demonstrate that potentials which 
match the observations within the central 200 pc of the Galaxy are crucial for understanding 
HVSs. Our approximations to the potential may be useful for other astrophysical problems 
which connect the central regions of the Galaxy to its outer reaches. 

We show that important aspects of the median speeds and shapes of the observable 
velocity distributions of HVSs are set by the Milky Way potential at r < 200 pc. For the 
potentials we construct to match the observations in this central region, median speeds 
are larger and velocity histograms are broader. Thus, low velocity ejected stars have less 
penetration into the outer halo at every stellar mass. These results indicate that HVSs might 
be useful probes of the triaxial potential of the central Galaxy. 

The models predict the spatial and velocity distributions of observable HVSs. They also 
provide a physical understanding of the origin of the dependences of these distributions on 
stellar mass for m = 0.6-4 M and on distance in the halo for r = 10-200 kpc. Here, we 
enumerate the main predictions of the model. We concentrate on the predictions resulting 
from our model potentials which match observations in the Galactic Center. For all of these 
issues, we consider the subtle effects on unequal mass binaries in the text (§3); here, we focus 
on the results for equal mass binaries. 

• Stellar evolution affects the observability of HVSs. It removes them at large radii where 
the travel time from the Galactic Center exceeds the stellar lifetime. 

Because the lifetimes of low mass stars are long, the observable properties of low mass 
HVSs depend only on the potential and the properties of binaries at the Galactic 
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Center. 

• For 3-4 M Q ejected stars, the fraction of unbound HVSs increases dramatically with 
r in the Galactic halo (Fig. [5]). Because the travel time from the Galactic Center is 
a significant fraction of the main sequence lifetime for these stars, the median stellar 
age of 3-4 M HVSs increases monotonically with r. 

In contrast, bound stars dominate the population of long-lived, low mass HVSs. Thus, 
the median stellar ages are essentially independent of r. 

• The velocity distribution of 3-4 M HVSs is asymmetric with a long tail toward positive 
velocities. The shape of the velocity distribution again reflects coincidence of stellar 
lifetimes and travel times. 

The velocity distributions of HVSs are increasingly symmetric with decreasing stellar 
mass. 

• As emphasized by Bromley et al. (2006), the short lifetimes of 3-4 M HVS require 
that they be injected with large velocities to reach the outer halo. Thus the median 
speed of massive HVSs increases with r. 

For low mass stars, the median speed depends only weakly on r. For all stars, the 
median speed increases with stellar mass at fixed r. 

• The density profiles of unbound HVSs are approximately p oc r~ n with n = 2-2.5. For 
the most massive stars, their finite lifetime removes stars at large r, steepening the 
density profile (n > 3 at r > 80 kpc). 

The density profile for bound, mostly low mass HVS is roughly p oc r~ 3 . We compute 
the detailed behavior of these profiles (Fig. H|). 

• Scaled density profiles show that the relative numbers of observable unbound HVSs 
as a function of stellar mass are relatively independent of the stellar lifetime. They 
depend mostly on the mass function at the Galactic Center at the time of ejection. 

In contrast, the relative numbers of bound HVSs are a function of the stellar lifetimes 
and the mass function. 

• We predict the relative observable numbers as a function of stellar mass for the range 
0.6-4 M . In a magnitude limited survey, the main factors that set detectability 
are (1) the accessible volume, (2) the fraction of unbound HVSs, (3) the asymmetry 
of the velocity distribution, and (4) contamination by indigenous stellar populations. 
Detection of HVSs is increasingly difficult with decreasing mass because all of these 
issues become less and less favorable. 
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We also argue that post main sequence stars are poor targets because their detectability 
is subject to yet another limit, the very short lifetimes of these phases. 

Samples of HVSs which are sufficiently large to explore these predictions provide a 
strong test of the model for injection of the stars into the Galactic potential. Coupled with 
observations of the stellar population at the Galactic Center, observations of HVSs provide 
promising probes of the bin ary populatio n, the stellar mass function, and the central potential 



of the Milky Way (see also IPeretsI 120071 ) 



So far, the observations by iBrown et al.l (l2007bl . and references therein) indicate that 



there are ~ 100 detectable HVS with m ~ 3-4 M . Dectection of comparable lower mass 
populations may be feasible. Samples of hundreds of HVSs promise strong tests of models 
like the one we construct. 

We thank J. Dubinski for helpful discussion on models for the Galactic potential and 
O. Gnedin and H. Perets for helpful comments on the manuscript. We acknowledge support 
from the NASA Astrophysics Theory Program through grant NAG5-13278. 
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Table 1. Adopted Properties of Main-Sequence Binary Stars 



Mass (M Q ) 


t ms (Gyr) 


Rrns (Rq) 


a min,q=l (AU) 


Q>min, 7712—0.6 (ATJ) 


0.6 


13 


0.6 


0.015 


0.015 


1 


10 


1.3 


0.032 


0.029 


1.5 


2.9 


2.5 


0.061 


0.051 


2 


1.2 


3.3 


0.081 


0.063 


3 


0.35 


4.7 


0.115 


0.084 


4 


0.16 


5.5 


0.134 


0.100 



Note. — These parameters are described in §2.1 of the text. 



Table 2. Parameters for Three-Component Galaxy Potentials 



Paper 


M b (10 9 M Q ) 


n (kpc) 


M h (10 12 M Q ) 


rh (kpc) 


M d (10 10 M Q ) 


a d (kpc) 


b d (kpc) 


DB98 


2.50 


0.05 


0.8 


21.8 


5.0 


2.4 


0.18 


GG05 


10.00 


0.60 


1.0 


20.0 


4.0 


5.0 


0.30 


YM07 


10.00 


0.50 


1.4 


25.9 


4.0 


5.0 


0.30 


here 


3.76 


0.10 


1.0 


20.0 


4.0 


2.0 


0.30 



Note. — These paramet ers are des cribed in §2.2 and E g. ITTI of t he text. The cited m odels are DB 
= Dehnen & Binnevl (Il998h . GG05 ^Cmedin etall (|200fih . YM07 = lYu fc Madaul (|2007h . and here = 
this paper 
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Table 3. Mass- Radius Relation for Milky Way 



Model 


M w (10 7 M Q ) 


M 120 (10 9 M Q ) 


M 210 (1O 9 M ) 


M 8000 (10 10 M Q ) 


M 1000 oo (1O 12 M ) 


Target 


3.0 


0.8 


2.7 


9.0 


0.7 


BKG06 


2.9 


1.2 


2.2 


9.0 


1.1 


here 


3.0 


1.1 


1.8 


8.2 


1.0 


DB98 


2.1 


0.4 


0.9 


9.0 


0.6 


WD05 


7.0 


1.3 


2.7 


9.0 


0.7 


GG05 


0.3 


0.3 


0.7 


8.2 


1.0 


YM07 


0.3 


0.3 


0.7 


7.8 


1.2 


SGM07 


2.3 


1.1 


2.4 







Note. — These parameters are descr ibed in §2.2 of the text. Referenc es are as in Ta bled except 
for W D05 = IWidrow fc Dubinskil (|2005l) . BKG06 = lBromlev et~ZI (|2006l) and SGM07 = IStolte et al 
(|2007h . 
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Table 4. Median Speeds of Ejected Stars at 30-60 kpc 



0.6 M 1 M 1.5 M Q 2 M q 3 M 4 M Q 



Simple spherical potential 



0.6 M 


7 


10 


32 


73 


211 


321 


1 Mq 


5 


7 


29 


71 


224 


351 


1.5 M 


5 


3 


32 


77 


229 


361 


2 M 


5 


5 


19 


77 


235 


364 


3 M 


1 


6 


15 


49 


238 


369 


4 M 


6 


3 


18 


39 


179 


372 


Three component potential from Gnedin et al. (2005) 


0.6 M 


5 


4 


21 


42 


135 


234 


1 M 


2 


4 


19 


44 


137 


238 


1.5 M 


-1 


1 


17 


42 


137 


240 


2 M 


3 


-1 


9 


43 


139 


241 


3 M 


1 


3 


11 


24 


135 


242 


4 M 





3 


9 


22 


104 


240 


Three component potential from this paper 


0.6 M 


7 


10 


30 


70 


199 


309 


1 M 


6 


7 


30 


70 


206 


330 


1.5 M 


4 


3 


31 


70 


209 


340 


2 M 


6 


7 


17 


69 


213 


342 


3 M 


1 


10 


16 


43 


216 


345 


4 M 


4 


4 


17 


39 


161 


347 



Note. - Velocity data (in units of km s _1 ) derived from 
simulations with the simple spherical Galactic potential (first 
set of entries ) , the three component Galactic potential from 
Gnedin et all |2005l ) (second set of entries), and the three com- 



ponent Galactic potential derived here (third set of entries). The 
bold entries correspond to the speed of the primary star in the 
binary progenitor; otherwise the speed corresponds to the sec- 
ondary star. 
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Table 5. Relative numbers of ejected stars 



Mass (M ) 


M v 


n e ,2i 


n e ,23 


n Ut 2i 


n u ,23 


0.6 


8.3 


160 


400 


8 


16 


1.0 


5.0 


125 


300 


6 


20 


1.5 


3.1 


55 


110 


5 


20 


2.0 


1.7 


30 


70 


8 


15 


3.0 


-0.3 


15 


15 


4 


4 


4.0 


-0.9 


1 


1 


1 


1 



Note. — These parameters are described in §5 of the main text. 
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Fig. 1. — Acceleration in the Galactic potential. Lower panel: decomposition of our Galaxy 
potential into contributions from the central black hole (thin solid black line), bulge (dotted 
line), disk (short dashed line), and halo (long dashed line). The thick solid line is the 
combined potential for comparison with the simple spherical po tential (thick grey line). 
Upper panel: as in the lower panel but for the iGnedin et al.l (120051 ) Galaxy potential. 
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Fig. 2. — Radial velo city profiles of ejecte d stars in different Galaxy potentials. The labels 
indicate profiles in the lGnedin et al.l (120051 ) potential (G), our three component potential (K), 
and the simple spherical potential (S). The differences in the three sets of profiles starting 
at Vq = 500, 750, and 1000 km s _1 illustrate the impact of the central potential in setting 
the observable properties of HVS in the Galactic halo. 



-40- 




Fig. 3. — Median age of ejected stars on their first pass through the halo as a function of 
Galactocentric distance. The labels indicate the masses of the ejected stars. 
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Fig. 4. — Relative density profiles of unbound stars ejected from equal mass binaries as a 
function of Galactocentric distance. In each panel, the labels for all ejected stars (solid lines) 
and for unbound stars (dashed lines) indicate the mass of the ejected stars. 
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Fig. 5. — Fraction of unbound stars as a function of Galactocentric distance for equal mass 
binaries. The mases are as in Fig. HI with 4 M for the top curve and 0.6 M Q for the bottom 
curve. 
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Fig. 6. — Predicted radial velocity distributions for equal mass binaries. The light grey and 
black histograms show the distributions in two different Galactocentric distance ranges; the 
clear histograms show results for all objects at 10-60 kpc. 
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Fig. 7. — Median speed as a function of Galactocentric distance for stars ejected from equal 
mass binaries. The labels indicate the stellar mass; the shade of the points increases with 
mass. 
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Fig. 8. — As in Fig. [7] but for outgoing stars. 
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Fig. 9. — As in Fig. [7J but for secondary stars in unequal mass binaries with a 4 M 
primary. 
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Fig. 10. — As in Fig. [7J but for 4 M Q primary stars in unequal mass binaries. 
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Fig. 11. — As in Fig. El but for stars in the three component and spheric al Galaxy pot e ntials . 
Upper two panels: velocity histograms in the plane of the disk in the I Gnedin et al.l (120051 ) 
potential for 0.6 M Q stars (left panel) and for 4 M stars. Middle two panels: As in the top 
panels for ejections along the z-axis of our three component potential. Lower two panels: As 
in the top panels for the simple spherical potential of Figure El 



10 



100 10 

r (kpc) 



100 



Fig. 12. — Relative density profiles for HVS as a function of stellar mass. In all panels, we 
normalize the density profiles for 0.6 M stars to unity at all r. Upper panels: profiles for all 
ejected stars for mass functions with q = -1.35 (left panel) and q = -1 (right panel). Lower 
panels: as in the upper panels for unbound stars. Compared to a population 0.6-1 M ejected 
stars, the relative abundance of more massive ejected stars declines monotonically with 
stellar mass (upper panels). For unbound stars (lower panels), 3-4 M Q (1.5-3 M ) stars are 
relatively more abundant than 0.6-1 M stars at r < 50-80 AU (r > 80-100 AU). Although 
these results show that 2-4 M ejected stars are factors of 10-100 less numerous than 0.6-1 
M ejected stars, magnitude limited surveys for HVS are likely to yield comparable numbers 
of massive ejected stars (see Table |5]). 



